Method of matched asymptotic expansions

In mathematics, particularly in solving singularly perturbed differential equations, the method of matched asymptotic expansions is a common approach to finding an accurate approximation to a problem's solution.

1 Method overview
2 Simple example
3 See also
4 References

Method overview

In a large class of singularly perturbed problems, the domain may be divided into two subdomains. On one of these, the solution is accurately approximated by an asymptotic series found by treating the problem as a regular perturbation. The other subdomain consists of one or more small areas in which that approximation is inaccurate, generally because the perturbation terms in the problem are not negligible there. These areas are referred to as transition layers, or boundary or interior layers depending on whether they occur at the domain boundary (as is the usual case in applications) or inside the domain.

An approximation in the form of an asymptotic series is obtained in the transition layer(s) by treating that part of the domain as a separate perturbation problem. This approximation is called the "inner solution," and the other is the "outer solution," named for their relationship to the transition layer(s). The outer and inner solutions are then combined through a process called "matching" in such a way that an approximate solution for the whole domain is obtained.^[1]^[2]^[3]

Simple example

Consider the equation

$\epsilon y'' %2B (1%2B\epsilon) y' %2B y = 0,$

where $y$ is a function of $t$ , $y(0)=0$ , $y(1)=1$ , and $0<\epsilon\ll 1$ .

Outer and inner solutions

Since $\epsilon$ is very small, our first approach is to find the solution to the problem

$y'%2By=0,\,$

which is

$y=Ae^{-t}\,$

for some constant $A$ . Applying the boundary condition $y(0) = 0$ , we would have $A=0$ ; applying the boundary condition $y(1) = 1$ , we would have $A=e$ . At least one of the boundary conditions cannot be satisfied. From this we infer that there must be a boundary layer at one of the endpoints of the domain.

Suppose the boundary layer is at $t=0$ . If we make the rescaling $\tau = t/\epsilon$ , the problem becomes

$\frac{1}{\epsilon }y''(\tau ) %2B \left( {1 %2B \epsilon } \right)\frac{1}{\epsilon }y'(\tau ) %2B y(\tau ) = 0,\,$

which, after multiplying by $\epsilon$ and taking $\epsilon = 0$ , is

$y'' %2B y' = 0\,$

with solution

$y=B-Ce^{-\tau}\,$

for some constants $B$ and $C$ . Since $y(0)=0$ , we have $C=B$ , so the inner solution is

$y_I = B\left( {1 - e^{ - \tau } } \right)= B\left( {1 - e^{ - t/\epsilon } } \right).\,$

Matching

Remember that we have assumed the outer solution to be

$y_O=e^{1-t}.\,$

The idea of matching is for the inner and outer solutions to agree at some value of $t$ near the boundary layer as $\epsilon$ decreases. For example, if we fix $t=\sqrt{\epsilon}$ , we have the matching condition

$\lim_{\epsilon \to 0} y_I \left( {\sqrt \epsilon } \right) = \lim_{\epsilon \to 0} y_O \left( {\sqrt \epsilon } \right),\,$

so $B=e$ . Note that instead of $t=\sqrt{\epsilon}$ , we could have chosen any other power law $t=\epsilon^k$ with $0<k<1$ . To obtain our final, matched solution, valid on the whole domain, one popular method is the uniform method. In this method, we add the inner and outer approximations and subtract their overlapping value, $\,y_{overlap}$ . In the boundary layer, we expect the outer solution to be approximate to the overlap, $y_O\,\sim\, y_{overlap}$ . Far from the boundary layer, the inner solution should approximate it, $y_I\,\sim\, y_{overlap}$ . Hence, we want to eliminate this value from the final solution. In our example, $y_{overlap}\,\sim\, e$ . Therefore, the final solution is,

$y(t) = y_I %2B y_O - e = e\left( {1 - e^{ - t/\epsilon } } \right) %2B e^{1 - t} - e = e\left( {e^{ - t} - e^{ - t/\epsilon } } \right).\,$

Accuracy

Substituting the matched solution in the problem's differential equation yields

$\epsilon y'' %2B \left( {1 %2B \epsilon } \right)y' %2B y = 0,\,$

which implies, due to the uniqueness of the solution, that the matched asymptotic solution is identical to the exact solution up to a constant multiple, as it satisfies the original differential equation. This is not necessarily always the case, any remaining terms should go to zero uniformly as $\epsilon \rightarrow 0$ . As to the boundary conditions, $y(0)=0$ and $y(1) = 1 - e^{1 - 1/\epsilon }$ , which quickly converges to the value given in the problem.

Not only does our solution approximately solve the problem at hand; it closely approximates the problem's exact solution. It happens that this particular problem is easily found to have exact solution

$y(t) = \frac{{e^{ - t} - e^{ - t/\varepsilon } }}{{e^{ - 1} - e^{ - 1/\varepsilon } }},\,$

which, as previously noted, has the same form as the approximate solution bar the multiplying constant. Note also that the approximate solution is the first term in a binomial expansion of the exact solution in powers of $y(1) = e^{1 - 1/\epsilon }$ .

Location of boundary layer

Conveniently, we can see that the boundary layer, where $y'$ and $y''$ are large, is near $t=0$ , as we supposed earlier. If we had supposed it to be at the other endpoint and proceeded by making the rescaling $\tau = (1 - t)/\epsilon$ , we would have found it impossible to satisfy the resulting matching condition. For many problems, this kind of trial and error is the only way to determine the true location of the boundary layer.^[1]

References

^ ^a ^b Verhulst, F. (2005). Methods and Applications of Singular Perturbations: Boundary Layers and Multiple Timescale Dynamics. Springer. ISBN 0-387-22966-3.
^ Nayfeh, A. H. (2000). Perturbation Methods. Wiley Classics Library. Wiley-Interscience. ISBN 978-0-471-39917-9.
^ Kevorkian, J.; Cole, J. D. (1996). Multiple scale and singular perturbation methods. Springer. ISBN 0-387-94202-5.